^ . Multiple scattering approach to elastic low-energy electron 
§ collisions with the water dimer 

! D. Bouchiha, L. G. Caron^ J. D. Gorfinkiel* and L. Sanche 

! Groupe en Science des Radiations, Departement de Medecine Nucleaire et Radiobiologie, 

[~^. ' Faculte de Medecine, Universite de Sherbrooke, Sherbrooke, Qc, J1H 5N4 Canada 

*Department of Physics and Astronomy, The Open University, Walton Hall, MK7 6AA 
Milton Keynes, UK 

<~| ■ ^Departement de Physique et Regroupement Quebecois sur les Materiaux de Pointe, 

Universite de Sherbrooke, Sherbrooke, Qc, J1K 2R1 Canada 

D ' Abstract. Multiple scattering theory is applied to low-energy electron collisions with a 



43 
Oh 



O 



complex target formed of two molecular scatterers. The total T-matrix is expressed in terms 
of the T-matrix for each isolated molecule. We apply the approach to elastic electron-(H20)2 



CO 

ZJ ' collisions. Following the method developed in our previous work on crystalline ice (Caron 

' 5/5 . et al. 2007), we impose a cut-off on the dipole outside the R-matrix sphere and an energy 

^"V dependent cut-off on the angular momentum components of the monomer T-matrix. An R- 

matrix calculation of electron-dimer collisions is performed in order to evaluate the accuracy 
of the multiple scattering approach. The agreement between the two calculations is very good. 
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1. Introduction 



t^- \ In biological systems, water is present in a wide range of environments. As a highly 
polar molecule, water tends to trap low-energy electrons (LEEs) and is expected to play an 
important role in LEE-induced processes (Garrett et al. 2005, Ptasinska & Sanche 2007). 
Electron collisions with water vapour being the simplest to study, extensive experimental 
(Lozier 1930, Schulz 1960, Compton & Christophorou 1967, Sanche & Schulz 1972, Belie 
et al. 1981, Curtis & Walker 1992, Fedor et al. 2006, Rawat et al. 2007) and theoretical (Gil 
et al. 1994, Morgan 1998, Gorfinkiel et al. 2002, Haxton et al. 2007) work has been carried out 
(for a comprehensive review of work on electron scattering from gas phase water up to 2004 
see Itikawa & Mason 2005). In the more complicated condensed phase systems, experimental 
studies of electron scattering from amorphous solid water, porous amorphous solid water 
and crystalline ice have been reported (Rowntree et al. 1991, Simpson et al. 1997, Michaud 
et al. 2003, Herring-Captain et al. 2005). In water clusters, formation of (H 2 0)~ by attachment 
of slow electrons has also been investigated (Knapp et al. 1986, Knapp et al. 1987, Weber 
et al. 1999, Barnett et al. 1989, Lee et al. 1991) in order to elucidate the mechanism of 
electron solvation. To our knowledge, no theoretical investigation of electron scattering with 
clusters or condensed water have been reported; available methods are mainly developed for 
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electron collisions with gas-phase molecules. It is therefore highly desirable to develop a new 
theoretical approach for condensed-matter and cluster environments especially since many 
electron-driven processes occur in these media. 

The main idea behind the multiple scattering (MS) method is to separate the potential 
of a complex target into non-overlapping regions with each region taken as a single scatterer. 
The impinging wave on each scatterer is composed of the incident plane wave and the wave 
scattered from the other scattering centres. The advantage of such an approach is that one 
can determine the cross section for large systems by combining information from its subunits. 
This enables the study of larger targets (macromolecules, molecular clusters, etc.) than is 
possible using the methods currently available. 

The MS approach to LEE scattering was previously proposed by two of the authors 
(L. Caron and L. Sanche) (Caron & Sanche 2003, Caron & Sanche 2004, Caron & 
Sanche 2005) and its use in conjunction with R-matrix calculations was recently validated 
(Caron et al. 2007) (hereafter referred to as I). The proposed approach combines scattering 
gas-phase data, obtained from accurate ab initio R-matrix calculations, to derive scattering 
information for a condensed-matter situation. The partitioning of the space in the R-matrix 
method makes it ideal for a multiple scattering approach of the "muffin-tin" type. The method 
is forcibly less rigorous than the ones available for small gas-phase targets, but it is fast and, 
for the systems tested so far, yields very good results. 

In I, we have used the MS theory in the condensed phase, using the equivalent KKR 
(Korringa, Kohn, and Rostoker) approach, in combination with gas-phase R-matrix data in 
the form of T-matrices. We were able to derive the band structure of a molecular crystal (ice). 
The present work constitutes a first attempt to determine elastic cross section (CS) for LEE 
scattering from a molecular cluster using the MS approach. The water dimer is an excellent 
choice for such an attempt, primarily because, making use of our previous experience in 
electron-H 2 collisions (Gorfinkiel et al. 2002), we were able to perform accurate R-matrix 
calculations of electron-dimer collisions and test the MS technique. In addition, water is 
present in a large variety of environments and information on LEE interactions with it is 
highly relevant. For example, the characteristics of the water in some of the solvation layers 
surrounding the double-helix structure of DNA are not that of bulk water (Becker et al. 1997), 
so may be explained in terms of clustered H 2 0. 

The main question we address in this paper is whether a multiple scattering approach is a 
good alternative to the computationally heavy ab initio methods for large to very large targets. 
For this purpose, we investigate how to apply the cut-offs found in I in order to perform a MS 
calculation of LEE collision with the water dimer. 

2. Theory 

In this section we present a multiple scattering method to calculate the total T-matrix, 
T tot , describing LEE collisions from a target consisting of two identical scatterers from 
their individual T-matrices T. The monomers, labelled with the index n, are located at 
R n = -R-i- = ii?. 
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The asymptotic wave function far away from the centre of mass of the two molecules can 
be expressed as follows: 

1 
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where r is the relative coordinate of the scattered electron and the centre of mass of the target. 
The Yl are spherical harmonics with L = (Z,m), ji denotes the spherical Bessel functions 
and h\ the spherical Hankel functions of the first kind, k = \[E with E the energy of the 
incoming electron and /° is an amplitude factor for the incident plane wave. The first term on 
the right side of (1) corresponds to the total incident wave and the second to the total scattered 
wave. 

We need to re-expand the incident wave in (1) around the centre of mass of each scatterer 
n. We use the general re-expansion formula from Dill & Dehmer (1974) which, for the 
spherical harmonics defined by Messiah (Messiah 1962), is 

Y L ,(r)h\t\Kr) = £ i h+h - l \-l) m ' F l ^ 2 _ m ,Y Ll (f n ) Y L2 (R n ) h&{ K r n )j b {kIQ (2) 
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where f n = f — R n and 
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is the Wigner 3-j symbol. 



It follows that the incident wave on each scatterer is: 

*f>i(r n ) = Y L'{r n )j v (Kr n )gl L , (4) 
v 

where 

g° nLl =Y. MV l lL fl (5) 

L 

or, in matrix form 

g° = Ml ■ f° (6) 

with 

Ml n LlL = Y,^ ll+h ~ l (-^ m F^X-mYL 2 (Rn)j k (KR n ) . (7) 

L 2 

We can now incorporate the multiple scattering between the two scatterers to obtain a 
global expression for the total impinging wave on each molecule. This has the form 

i>i(r n ) = Y L {r n )ji(Kr n )g nL . (8) 

L 

One can identify our g nL with the expression Ani l e l k ■ R n B^ in the analysis of Caron & 
Sanche (2004). By comparing Eq. (8) to the total impinging part of their Eq. (1), i.e. the first 
term on the right-hand side, one then gets from their Eq. (2) 
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with i? nn / = R n — R n i, n' = —n. 

X7l> = \ E^ + '^(-l) m ' 2 <5,- m ^L a (£nn') Mn.') (10) 

expresses how a scatterer influences its neighbour. Note that for n = n' Xll> = 0- The second 
term on the right side of (9) is the total scattered wave from the second scatterer n'. In matrix 
form, (9) can be written as 

g=(I- X T)" 1 -g°. (11) 

Now, we have to combine the amplitudes scattered from each scatterer 

^(r) = E Mr n ) = E E Y L ir n )hf\Kr n )T2, L g nL . (12) 

n n LL' 

around the common origin (centre of mass) using 

Y v {f n )hf\nr n ) = E i h+l2 ~ l '(-ir'^'!i -m'^ (r) Y L2 (-R n ) h^\nr) 3l2 (nR n ) .(13) 
One gets 

M^=J2Y Ll (r)h^(Kr) E M2 Li l,T2, L g nL (14) 

Li A n,L,L' 

where 

M2 Ll n L , = E i h+l2 - l \-^) m ' F^'X^ m> Y L2 (-R n ) j l2 {kRJ . (15) 

Combining (14), (6), (11) and referring to the general expression (1), we have 

T tot = M2-T-(I-xT)~ 1 -Ml . (16) 

We thus obtain the T-matrix for two identical scatterers from their individual T-matrices. 
Ml and M2 account for the transformation from the monomer to the dimer' s centre of mass 
and (I — %T) _1 expresses the multiple scattering between the two scatterers. 



3. Characteristics of the calculations 



No experimental or theoretical data on elastic LEE collisions with the water dimer are 
available. Therefore, in order to test the quality of our MS results, we have performed an 
ab initio scattering calculation on the gas-phase dimer using the R-matrix method and the UK 
polyatomic R-matrix suite (Morgan et al. 1998). 

The gas-phase water dimer in its ground state equilibrium geometry has C s symmetry 
(Figure 1). The geometry parameters from Park C-Y et al. (2001) were used in the R- 
matrix calculation. The geometry of the water monomer is very slightly changed upon dimer 
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formation (see Table 1). In the MS calculation, the dimer was formed by putting together two 
water molecules separated by d ~ 5.5 bohr along the Z-axis. This means that there are very 
minor differences between the geometries used in the R-matrix and MS calculations. The 
most notable is the hydrogen bond r OH which is at most ~ 0.013 bohr longer in the dimer 
than in the monomers we use to build the MS results. 

3.1. R-matrix calculation for the dimer 

A detailed description of the R-matrix method as applied to polyatomic molecules within 
the fixed-nuclei (FN) approximation can be found in Morgan et al. (1997) and Morgan et al. 
(1998). The method is based on splitting coordinate space into two regions separated by a 
sphere, henceforth called R-sphere, centred on the centre of mass of the molecule. The radius 
of the R-sphere is chosen in such way that the electronic density of the target is negligible 
outside it. As a consequence, exchange and correlation effects can be neglected in this outer 
region and a long-range multipole expansion used to represent the electron-target interaction. 
Inside the R-sphere, both effects are significant and are therefore taken into account using 
rigorous quantum chemistry methods. The wavefunction for the target + electron system is 
then expressed in terms of a close-coupling expansion. 

The calculation was performed following previous work on electron scattering from 
isolated H 2 (Gorfinkiel et al. 2002). We used the equilibrium geometry from Park C-Y et al. 
(2001) and the same basis set (including the diffuse functions) employed for the monomer 
work by Gorfinkiel et al. (2002). Since in this work we are concerned with the elastic 
scattering, only the ground electronic state was considered and included in the close-coupling 
expansion. We therefore produced natural orbitals exclusively from this state and then used 
them in a CASCI (complete active space configuration interaction) calculation in which 8 
electrons are frozen; this generated around 7000 configurations. With this model, we obtained 
a good value for the ground state energy: -152.18 hartree, compared to -152.67 hartree from 
the most accurate calculation. We also obtained excellent agreement with the experimental 
dipole moment: 1.065 a.u. in our calculations compared to 1.041 a.u.. 

We used an R-sphere radius of a=13 bohr. In order to confirm that all the electronic 
density was contained inside the sphere, tests were performed for radii of a=14 and a=15 bohr. 
Only small differences at very low energies were found between the cross sections calculated 
with the three radii (Gorfinkiel 2008). The continuum orbitals describing the scattered 
electron were expanded in a basis of GTOs with I < 4 centred on the centre of mass of 
the water molecule. 

3.2. Multiple scattering calculation 

In I, we showed that gas-phase R-matrix data can be efficiently used in conjunction with MS 
theory in a non-gaseous environment by deriving the band structure of a molecular crystal 
(ice). The fundamental lesson learnt from this KKR study is that for the multiple scattering 
term, we need a trimmed T matrix T c , without dipole contribution for radii r > a c and with 
an angular momentum cut-off. Thus two different monomer T-matrices are needed: T c for 



Multiple scattering approach to low-energy electron collisions with the water dimer 



6 



the multiple scattering term (I — %T C ) and a matrix T dip for which no cut-offs are applied. 
With this in mind Eq. 16 can be rewritten: 

T = M2 T dip ■ (I - xTc)- 1 ■ Ml . (17) 

The first step in our multiple scattering calculation involves generating the monomer T- 
matrices T c and T d;p required in equation (17). The description of LEE scattering by an 
isolated H 2 molecule is the same than in I. Since a detailed description of the calculation 
was given there, we will limit ourselves here to a brief summary. We have used the R-matrix 
method within the fixed-nuclei approximation. For the target description we used the DZP 
basis set (Dunning 1970) for the Oxygen and a TZ basis (Dunning 1971) for the Hydrogen. 
A CASCI model was used with only 2 frozen electrons and the calculation included only the 
ground state. The scattering calculation was performed with an R-sphere radius of a = 6 bohr. 
A GTO basis for the continuum with / < 4 was generated for this radius using the program 
GTOBAS (Faure et al. 2002) . The obtained T-matrices were rotated from the R-matrix 
coordinates to their corresponding position in the dimer. The procedure to do so was explained 
previously (Equation (24) in I) and the reader is referred to that work for further details. 

As mentioned above, two cut-offs are applied when generating T c . A first cut-off on 
the range of the molecular dipole is achieved by removing the dipole field for r > a c . 
The cut-off radius a c depends only on the dipole moment of the molecular monomer and 
should therefore be the same for the ice and the dimer calculations. In I, we have found that 
a c ~ 6.5 bohr yields the best value for the electron's effective mass in ice. An analysis of 
the electron-(isolated) water cross section obtained for different values of a c supports this 
finding (the reader is referred to Fig. 1 in I): for too small cut-off radii (a c < 6), the high 
energy cross section is unphysical whereas for higher radii (a c > 6) the cross section exhibits 
a dipole-driven behaviour (it increases rapidly at low energy). Although this behaviour is 
physical and should be expected, it is the aim of the cut-off to at least partially eliminate it. 
This can be understood as an attempt to ensure that T c involves only the contribution of the 
target potential in the small region around the scatterer. That leaves a c ~ 6 bohr as the most 
satisfactory radius. 

The second cut-off, in the angular momentum components I of the scattering matrix, is 
critical and depends both on the electron kinetic energy E and the intermolecular distance d. 
The restriction is related to the angular momentum energy barrier E(l, r) = 1(1 + l)/r 2 . The 
KKR calculation predicts that only angular momentum values I < l c should be retained such 
that E(l c , d) < E e < E(l c + 1, d) with E(l, d) = 1(1 + l)/d 2 Ry. That is, only electrons that 
scatter from one monomer with energy larger than E(l,d) can reach the other one and undergo 
multiple scattering. In the dimer, the two H 2 molecules are separated by d ~ 5.5 bohr 
(almost the same intermolecular distance than in crystalline ice in I). We expect then to have 
the same cut-off criteria than in the KKR calculation. 

Once the T-matrices for the monomer are calculated, the dimer T-matrix in equation 17 
is built. From it, we build the R-matrix at r = adi P = 9 bohr using the approach given in 
Appendix A and propagate it outwards (Baluja et al. 1982) using the exact dipolar field of the 
dimer. This is necessary in order to incorporate the effect of the dipole moment of the dimer 



Multiple scattering approach to low-energy electron collisions with the water dimer 



7 



in the MS description. To choose a dip we propagate the R-matrix of the dimer (obtained with 
the R-matrix calculation) inwards from the asymptotic region and then calculate the CS from 
it. The low-energy part of this CS decreases the more of the dipole field we remove (the 
further in we go) and takes an unsuspected and incorrect upturn for a radial distance of 9 bohr. 
We believe that for this distance the multipolar field of the dimer is no longer describable 
by its total dipole field. We finally generate the final T-matrix T tot which takes the dimer's 
dipole contribution into account. This T-matrix is a function of energy and of the angular 
momentum cut-off l c used for T c in the multiple scattering term of Eq. 17. The integral 
elastic cross section is then obtained using the well known formula: 

K IV 

In Figure 2, we present the CS for different l c values. As mentioned earlier, the partial 
wave expansion must be limited as the collision energy decreases. For example, for electron 
energies between 5.4 < E < 9 eV, the KKR calculation predicts that we can only include 
/ < 3 in the MS calculation. It can be easily seen from Figure 2, that the CS is ill behaved 
in this region (below ~ 7 eV) when l c = 4. A closer look at the CS for the other values of 
l c shows the presence of a critical energy below which the CS does not exhibit the right 
behaviour. These critical energies are at slightly lower energies than the value given by 
E(l c ,d). 

If we scrupulously restrict the calculation to integer values of l c for electron energies 
between E(l c ,d) and E(l c + l,d), the CS will be jagged i.e. discontinuous at E(l c ,d). 
Smoothing is needed. A two-point interpolation on the parameter I between l c and l c + 1, 
where I is the solution of E = 1(1 + 1)/ d 2 , would get rid of the discontinuities. But this would 
admix the CS at l c with any singular part of the CS at l c + 1. A quick look at Figure 2 indicates, 
for instance, that at 6 eV when l c = 3 the interpolation would spuriously bring in a sizable 
part of the huge peak of the l c = 4 CS. The procedure can be regularised by upward shifting 
E(l c , d) using a parameter 7 such that 7-E , s (/ c , d) = l c (l c + l)/d 2 . A value of 7 = 0.75 would 
shift the cut-off energies to the new values shown in Figure 2, beyond the threshold of the 
singular behaviour of the CS for l c + 1. The interpolation procedure can now be done safely 
with I a solution of ^E s = 1(1 + l)/d 2 . The CS is linearly interpolated between CS values 
corresponding to the closest lower and larger integer values of I. 

4. Results 

The elastic electron-(H 2 0)2 cross section versus impact energy is plotted in Figure 3. As a 
test of the accuracy of the MS calculation, we compare our results to the CS calculated with 
the R-matrix codes. The agreement between the two calculation is very good. The multiple 
scattering cross section is slightly higher than the R-matrix one but remains within 5% of 
the the fully ab initio CS for energies larger than 2.5 eV. We believe these differences are 
well within the range of usual experimental errors and ab initio calculation uncertainties. For 
targets with large dipole moments, it is customary to add to the R-matrix cross section a Born 
based correction (Chu & Dalgarno 1974) to account for the partial waves not included in the 
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ab initio calculation. This correction is the same for the R-matrix and MS calculations and 
for this reason has not been included in the cross sections plotted in Figure 3. 

Based on the results of I, the need for a "trimmed" T-matrix (T c ) in the MS term of 
(17) is obvious. Figure 2 illustrates the drastic impact of a larger angular momentum basis; 
including higher I at low energies introduces resonances and/or premature divergences in the 
CS caused by the resonant nature of the MS (I — xT) 1 term in Eq. 17. The cut-off in the 
long-range dipolar interaction is also very important. In I, we have shown that even though 
the value of the cut-off a c is not critical for band structure calculations, a too large or too 
small cut-off radius has a noticeable effect on the value of the effective mass (which was used 
for calibration purposes). The two cut-offs are necessary and have to be applied only to the 
MS term. The use of a full T-matrix T dip in (17) is, however, essential to get the right CS 
behaviour at low energies and to properly account for the potential produced by the dipole of 
a H 2 molecule. The need for Tdi p to include the dipole contribution arises from the fact 
that for r < a^ ip , the scattering electron feels the individual dipoles of the isolated water 
molecules. 

The radius a<a v has to be chosen in a proper way. It must allow the separation between 
the regions where the individual water dipoles dominate the electron-molecule interaction 
(r < ddip) and the region where the total dipole of the dimer dominates (r > a^p). Based on 
a geometrical analysis (Figure 4), one can see that a radius of 9 bohr is an adequate choice. In 
order to check the appropriateness of this value, we have tried several R-sphere radii ranging 
from 6 to 13 bohr and confirmed that 9 bohr yields the best cross section. 

It is evident from Figure 4, that a choice of a ~ a c ~ 6 bohr in our calculations does 
not exactly correspond to a non-overlapping potentials picture. For the spheres not to overlap 
we would need a radius around d/2 ~ 2.7 bohr. However, using such small radius does not 
provide a good representation of the target. This was already clear in I, where we established 
that small radii introduce a discontinuity at the boundary of the sphere that is too important to 
lead to good results. 

5. Conclusions and Discussion 

We have developed a technique to derive the total T-matrix of two scatterers from their 
individual T-matrices as described by Eq. (16). The technique first involves finding the T- 
matrix describing the scattering for an isolated molecular target. These are then combined 
using Eq. 17. Our main conclusion is that two different T-matrices, with different constraints, 
are needed to represent the electron-dimer interaction: two cut-offs have to be applied only to 
the MS term (T c but not Tdi p ). From the total T-matrix, the R-matrix at a^p is determined 
and then propagated to an asymptotic distance using the dipole moment of the cluster. This 
propagation is fundamental in order to incorporate the electron interaction with the true dimer 
dipole moment. One can argue that the need for T dip and a^ ip is a surface effect not present 
in the infinite crystal situation in I. In fact, the dimer is mostly surface. 

Our results show that a multiple scattering calculation can efficiently replace a standard 
ab initio scattering calculation for this system. Very good agreement is found between the MS 
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and the R-matrix cross sections. The MS method should allow us to calculate elastic cross 
sections for small water clusters that are nevertheless too big to be studied ah initio. The same 
value of a c should be employed and a d i P could be determined using the same geometrical 
considerations. How large a cluster could one attempt to treat using this procedure? The 
cluster should have a large surface to volume ratio. The arguments presented in Appendix B 
indicate that the number of water molecules should be much less than 500. That leaves a lot 
of room for fair sized clusters. 

For very large clusters, having a low surface to volume ratio, the general picture would 
be as follows. The physical space should be divided into three parts. In the innermost part, 
the bulk of the cluster, MS theory using T c in a uniform background optical potential U op , 
measuring the average polarisation potential energy between the muffins, would be used as 
in I. Note that the resonant MS (I — xT c ) _1 term of Eq. (17), which has been freed of any 
spurious peaks through the angular momentum cut-off algorithm, will yet show peaks that 
correspond to damped electronic modes of the cluster. These are, for instance, the modes that 
were calculated for crystalline ice in I. Then there would be a narrow surface region which 
adapts the bulk value of the optical potential to the vacuum level and in which T c would be 
used for the MS part and Tdi p for the exit part. Finally, there would be an outer region in 
which the R-matrix of the two inner regions can be propagated to infinity using the multipolar 
potential of the cluster. 

We have thus far discussed only elastic scattering. How can one deal with inelastic 
collisions? Although this is unchartered territory, the scenario will likely go as follows. In 
such a situation, new energy channels open up. For each of these, there will be energy- 
diagonal T-matrices T d (Ei) and cross-energy ones T nd (Ei, Ej). These will combine into a 
super-matrix T(Ei, Ej) = Td{Ei)5E i ,E j + T nd (Ei, Ej). The rest would basically be a repeat 
of what we have done for a single channel. One would define the MS part T c with cut- 
offs on the range of the dipole and on the angular momentum basis. Note that the latter 
is different for each energy channel. Moreover, as the energy increases, larger and larger 
angular momenta will be needed. Some adjustments on a c might also be required in order 
to include inelastic channels involving more extended molecular wavef unctions. Formally, 
Eq. (17) still applies to these super matrices. The resonant MS (I — xT c ) 1 term would not 
only exhibit peaks related to the damped electronic modes of the cluster but also structures 
caused by the inelastic contributions to T c which will in turn modulate the inelastic peaks 
in T dip . Interpolation of the total CS would now involve linear interpolation of the type 
described above in a multi-dimensional energy-channel space. The computer time consuming 
aspect comes from repeated cluster calculations of Eq. (17) for many values of the incoming 
electron energy Ei and l c (Ei). 
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Appendix A. Relation between R and K matrices 



We suppose in the following, that the R-matrix defined at a radius r = a describes the inner 
region r < a and that the potential outside the sphere is zero. For elastic scattering we have 
(Burke & Berrington 1993), 



F L {a)=Y,R L u r—F L , . (A.l) 

L' L J r=a 

with F L the reduced radial wavefunctions. 

For r > a, the wavefunction representing the scattering electron is 

^(0 = E— (A.2) 

L r 

In the region r > a, one can write 

F L (r) = fiL jifar) + f 2L (A.3) 

where j(z) = z j(z) and fj(z) = z rj{z). Fl(t) is a solution of energy E = k 2 /2 and angular 
momentum L of the Schrodinger equation in a zero potential. Substituting (A.3) in (A.l), we 
have 

f\L + fiL Ve(f^a) = anJ^Rw [fw fei^a) + f 2L > Vi>( Ka )] , ( A - 4 ) 

v 

where j[{z) = dje(z)/dz and fj'i(z) = df\t{z)jdz. One can rearrange this equation and write 
it in matrix form 



(N - a^RN'j f 2 = - (J - auRJ'j f x (A.5) 

where Jll 1 = 3i{i^o)8ll', N ll > = fje(Ka)SLu and /i, f 2 are the amplitude column vectors. 
This gives us: 

f 2 = - (N - a K RN') _1 (j - okRJ') h (A.6) 
We introduce the K-matrix, 

K = (N - a^RN')" 1 (j - a«RJ') . (A.7) 
Substituting (A.7) in (A.6) and then in (A.3) and (A.2), we have for r > a 

^(r) = (nr)- 1 J2 fn []i(Kr)5 LL , - K UL f] v {nr)] Y L ,{r) (A.8) 

LL 1 

and thus 

lim ^(r) = {nr)- 1 ^ f 1L [sm(nr - £n/2)5 LU + K UL cos(fi:r - £V/2)] F L ,(f) ,(A.9) 

which is conform to the definition of the K matrix. 

Eq. (A.7) gives a 'passage' for the transformation R <-> K. It is worth noting that the 
R-matrix depends on the radius but not the K-matrix (nor the transmission matrix T, or the 
scattering matrix S). So once we have the R-matrix for any radius a we can build the K-matrix 
using Eq. (A.7) and inversely. The T-matrix is then obtained from the K-matrix using the well 
known relation: 

T = 2iK(l — iKy 1 (A.10) 
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Appendix B. Surface to volume ratio 

One can estimate a surface to volume ration by the ratio of the number of surface molecules N s 
to the number of bulk molecules N b . One can get an estimate of these numbers by considering 
a spherical cluster of radius R c . The surface molecules are located in the shell between R c 
and R c — d where d is, as above, the distance between molecules. If the mid-shell surface area 
S = <i7r(R c — d/2) 2 is occupied by molecules forming a close packed array, the area occupied 
by one molecule would be S mo i ~ d 2 a/3/2 and thus N s = S/S mo i- Assuming the bulk of 
volume V = Att{R c — d) 3 /3 is occupied by close packed molecules occupying a volume of 
Vb = 0.93d 3 a/3/2, one gets N v = V/V&< A high surface to volume ratio can thus be quantised 
by the quantity SV = N s /N b > I. With d = 5.5 a.u. one gets R c < 4.6. One thus deduces 
that the cluster must contain much less than N s + N v « 500 molecules to have a large surface 
to volume ratio. 
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Tables and table captions 





H 2 


(H 2 0) 2 


*OHl 


1.81 


1.810 


?OH2 


1.81 


1.823 


0i 


104.5 


104.5 


r OH3= r OH4 




1.814 






5.497 






104.6 


a 




4.7 


P 




55.1 



Table 1. H2O and (H2CO2 geometries used in our calculations. The values of the dimer 
parameters (see figure 1 for identification) are from Park C-Y et al. (2001). Length units are 
in bohr and angles in deg. 
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Figure captions 
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Figure 2. Trial elastic electron-water dimer cross sections calculated using different angular 
momentum cut-offs l c . The MS term in Eq. (17) introduces unphysical resonances and/or 
premature divergences in the CS. The filled arrowheads point to the cut-off energies E(l c , d) 
while the empty arrowheads point to the shifted regularised energies E s (l c ,d) for l c = 
1, 2, 3, 4 respectively, from left to right. 
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Figure 3. Elastic electron-(H20)2 cross sections. The MS cross section was obtained 
following the smoothing prescription explained in the text. Notice that he first vertical 
excitation threshold is located around 6.99 eV (Valenzano et al. 2005). 




Figure 4. Radii used in the MS calculations; solid line: a<n v = 9 bohr, for inclusion of the 
dimer's dipole effect. Dashed line: R-sphere with radius a = 6 bohr used in the R-matrix 
calculation for the monomer, n and n' define the positions of the centre of mass of each water 
molecule. 



